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Abstract 

We compute the pion quark Generalised Parton Distribution H and quark Double Distributions in 
a coupled Bethe-Salpeter and Dyson-Schwinger approach in terms of quark flavors or isospin states. 
We use analytic expressions inspired by the numerical resolution of Dyson-Schwinger and Bethe- 
Salpeter equations. We obtain an analytic expression for the pion Generalised Parton Distribution 
at a low scale. Our model fulfils most of the required symmetry properties and compares very 
well to experimental pion form factor or valence parton distribution function experimental data. 
In addition, we have highlighted limitations of the so-called impulse approximation, which breaks 
symmetries when computing the valence parton distribution function. Doing so, we introduced new 
terms which were neglected before. Finally, we also shed light on a specific property of the pion 
GPD: Polyakov soft pion theorem. 


1 Introduction 


Generalised Parton Distributions (GPDs) were introduced in the 1990 and since then have been 
deeply studied both theoretically (see e.g. the reviews and experimentally. New results from 

Jefferson Laboratory facilities (JLab) on Deep-Virtual Compton Scattering (DVCS) on a proton target 
have recently been presented [T], [^. Several phenomenological parametrisations have been developed 
in order to fit the proton GPDs [9l4l4l|. Yet until now, none of them has been fully derived from QCD 
dynamics only. 

In order to generate the parton structure of hadron in a dynamical way, one can turn to the 
Dyson-Schwinger (DS) and Bethe-Salpeter (BS) equations. Those equations are coupled and can be 
solved using specific truncation schemes. In the following, we focus on the pion GPD, which is simpler 
to compute in the present framework than the proton one. The pion GPD inspired many theoretical 
studies [iiH in spite of the restricted set of existing experimental data, related to the valence Parton 


Distribution Function (PDF) 21[ and the form factors 2^,12^. Our approach to the pion GPD is 


detailed in the first section. The second section deals with the restoration of the symmetry under the 
exchange x I — x in the forward case. 
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2 GPD Computations 

The pion quark GPD is formally defined as: 
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where f = and ^ = —■^+- In our approach [2^, [2^, we define the pion GPD by the set of its 
Mellin moments M„(^, t): 

= J dx (2) 

These moments are computed in the triangle diagram approximation illustrated on figured! Consequently, 




Figure 1: Mellin moments computed in the triangle approximation for the pion valence dressed quarks. 
Left panel: quark case. Right Panel: anti-quark case. 


the inserted operators, depicted on figured! by a cross, are the local twist-two quark operators arising 
from the operator product expansion. They can be computed as: 
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iy -n S{k + y) iT^{k + ^,P + ^) S{k - P), 
where the Bethe-Salpeter vertices T^ and propagator S are taken as: 

S{p) = [-ij ■p + M]Am{p‘^), 

= TW' 
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These functional forms have been introduced in Ref. and allow one to fit the numerical solutions 
of the Dyson-Schwinger and Bethe-Salpeter equations. Here we will fit a single free parameter, M, the 
second one, u, being fixed to 1. This value of u is the one giving back the pion asymptotic Distribution 
Amplitude (DA). From that point, the GPD can be analytically reconstruct from its Mellin using the 
so-called Double Distributions (DDs). Indeed, performing the computation, it is possible to identify the 
two DDs a, t) and G^{j5, a, t) which are related to the pion GPD through the Radon transform: 



l-x 


Figure 2: Left Panel: reconstruction of the GPD when t —>■ 0. Support and continuity properties are 
well illustrated on this plot. Right Panel: PDF large-x behaviour, which fits the one predicted in the 
pertubative case. 

to compute algebraically the GPD provides us a significant advantage to check the different intrinsic 
properties that the model should respect. Due to the property of the algebraic DDs identified from our 
calculations, it is possible to conclude that the support property of the valence GPD, the continuity 
on the lines x = and the parity in ^ are fulfilled by our GPD model. Moreover, the DD formalism 
ensure the polynomiality property by construction. Finally, one can check that the resulting PDF gets 
the same large-x behaviour than the one predicted by perturbation theory. 


3 Soft Pion theorem 


In Ref. 27|, it has been shown that, it was possible to relate the pion DA to the pion GPD. In the 
kinematic limit when t —>■ 0 and ^ 1, one actually gets: 
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(9) 


where (p-^ is the pion DA. The algebraic model coming from equations (jlD-dT]) have been developed in 
order to describe And indeed, when computed in Ref. [2^, it leads to the asymptotic pion DA. 
Therefore, one could expect that, due to the soft pion theorem, when f goes to 1 and t vanishes, the 
algebraic GPD model tend to the pion asymptotic DA. It is not the case due to a specific feature 
of the algebraic model. Indeed, working only with the building blocks of the solution of the DS-BS 
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equations leads to a violation of the Axial-Vector Ward-Takashi identity (AVWTI), which gives in the 
chiral limit: 

p^n^ik, P) = ^- l^s-\k - + il5S-Hk + |)) , (10) 

where P is the total momentum entering the axial-vector vertex, k the relative one, and the r* the 
Pauli matrices. When computed consistently in the rainbow ladder truncation scheme, the full solution 
of the DS-BS equation fulfils the AVWTI. Therefore, inserting the reduction formula: 

2UT^{k]p) pf"T5f,{k,p), ( 11 ) 

2r^r^(A:;p) ^oc° p‘^T5{k,p), (12) 

allows one to write the triangle diagram in terms of axial-vector and pseudo-scalar vertices. Then, as 
shown in Ref. [^, most of the gluon ladders coming from the insertion of the AVWTI compensate 
each others, leading to the pion DA. Consequently, in the rainbow ladder truncation scheme, the 
soft pion theorem is automatically fulfilled providing that the Dyson-Schwinger and Bethe-Salpeter 
equations are solved consistently with the AVWTI. 


4 Additionnal contributions 


If it is possible to recover the soft pion theorem through a triangle diagram computation, such an 
approximation still has its own limitations. Focusing on the forward limit {i.e. t = 0 and C = 0); the 
PDF can be written as: 

= — ((30 — 15x -k — 2x^)x^ logx -k (3 -k 2x^)(l — x)^ log(l — x) (13) 

25 V 

-k (3 -k 15x -k — 2x^)x(l — x)^ . 


If this PDF is in agreement with the pertubative prediction at large x, it suffers a significant drawback, 
as it is not symmetric with respect to the x -H- 1 — x exchange. The small asymmetry is due to the 
fact that contributions have been neglected, due to the triangle approximation. Indeed, it has been 
shown in Ref. 29|], that one has to add new terms in order to fulfil the x -f-)- 1 — x symmetry: 
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It is then possible to compute and additional component to the PDF q^{x): 

72 / 

+ 4x + 9 ) (x — 1)^ log(l — x) + x^(2x((x — 3)x + 5) — 15) log(x) 
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(16) 


Then summing the different contributions, one gets the total PDF ql^\x): 
ql°\x) = q^{x)+qf{x) 


72 


x^(x(2x — 5) + 15) log(x) — (2x^ + x + 12) (x — 1)^ log(l — x) 


—2x(x — l)((x — l)x + 6) , 


(17) 


which fulfils the x 1 — x s ym metry as shown on figure ElThis analysis has also led to a new GPD 


Ansatz described in Ref. [28|, l30|, 



X 


Figure 3: Contributions to the PDF and their sum. 


5 Conclusion 

We haye deyeloped a new GPD model. It is based on algebraic building blocks of solutions of the 
Dyson-Schwinger and Bethe-Salpeter equations. Within this procedure, our model fulfils most of 
the required theoretical constraints. In addition, our model is in good agreement with available 
experimental data. Computations has been done analytically up to the end, allowing one to identify 
limitations of the triangle diagram approximation, and thus to improve it. We have also emphasised 
the key-role of the AVWTI in the realisation of the soft pion theorem. Future models of pion GPDs 
will therefore have to deal with those additional constraints. 
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